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FIELD EMISSION CHARGING OF METALLIC COLLOIDS 


By N. Stankiewicz 


SUMMARY 

With a free-electron model of metallic colloids, the charging time in an 
applied electric field is calculated for the alkali-group elements. A theoreti- 
cal thin-wire ionizer is discussed with regard to required operating potentials, 
necessary geometrical dimensions, and current -density capabilities. 


INTRODUCTION 

The theoretical improvement of electrostatic-thrust or efficiency resulting 
from the use of charged particles of large mass has been recognized for some time 
(ref. l). Current studies on heavy-particle thrustors have developed in two 
areas. One is concerned with employing large molecules and the other with ag- 
gregates of atoms or molecules of colloidal size (simply called colloids). Both 
types of heavy particles present the problem of a mass distribution in the beam 
resulting in a loss of thrustor efficiency. For colloids, however, reference 2 
indicates that the mass distribution can be quite narrow. 

At present, several methods of charging colloidal particles are being in- 
vestigated. The method to be considered in this report is that of ionization by 
field emission of an electron. For reasons that are elaborated in the next 
section, only metallic aggregates, and, in particular, only those of the alkali 
group, are discussed in detail. 

The field ionization method was studied in the investigation reported in 
reference 3 for the case of a gas of alkali atoms. Results were not encourag- 
ing, however, because of the excessive voltages required. This study was under- 
taken to determine if requirements might be less stringent for colloids than for 
atoms. 


GENERAL DISCUSSION OF COLLOIDS 

The structure of a colloidal particle is not clearly under stood. It is 
highly improbable that any well-ordered arrangement of atoms prevails such as in 
a crystal. It is more likely that the colloid resembles a liquid droplet, within 
which the atoms are randomly arrayed. 

The electronic structure of a colloid is similarly not understood. With the 



assumption that the density of a colloidal particle is equal to its bulk density, 
it would be expected that metallic colloids contain unbound electrons; that is, 
electrons shared by all the member atoms of the aggregate. This follows from 
the nearness of neighboring atoms in a solid or liquid and also from the low 
ionization potentials of metal atoms. On the other hand, atoms with high ion- 
ization potentials, in the solid or liquid state, would be devoid of free elec- 
trons and, therefore, would be nonconducting. Removing electrons from such 
aggregates would hence be equivalent to ionizing the constituent atoms and would 
be a much more difficult task because of the higher ionization potentials in- 
volved. 

To justify the hypothesis of free electrons, it may be pointed out that even 
a liquid metal is a good conductor despite the disorder of its atoms. An alter- 
native statement of the f ree-electron condition is the fact that the energy 
levels of the valence electrons are not determined solely by the atoms in which 
they were originally bound, but by all the atoms in the colloid; that is, the 
allowed energy states must be determined by solving for the eigenvalues of the 
wave equation in which the electrons are shared by all the ionized atoms in the 
aggregate. 

Because of the great difference in mass between electrons and atoms, the 
atoms can be considered as quasi -stationary, while the electrons move among them. 
Colloids, however, have no periodic structure, and thus an already difficult 
calculation for the energy levels becomes unmanageable. 

With the assumption that these allowed energy states E in the droplet are 
somehow determined, Fermi statistics can be formally applied to obtain the fol- 
lowing distribution function for the electrons: 


f(E) 
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(Symbols are defined in appendix A. ) 


The Fermi energy E 0 , however, would be both space and time dependent be- 
cause of the random spacing of the metallic ions and the slowly changing inter- 
ionic distances. These distances are slowly changing with respect to the move- 
ments of the electrons but are rapidly varying with respect to the motions exe- 
cuted by the colloid as a whole. Because of this, it would be expected that an 
average Fermi level could be assigned in equation (l). In lieu of detailed 
analysis and since this paper intends to give results indicative of the charging 
mechanism, the Fermi level is taken as equal to the Fermi level found in bulk 
materials. Elements with the smallest work functions will therefore be the 
easiest from which to remove electrons and thereby produce positively charged 
colloids. In view of this, the discussion is restricted to the alkali metals. 
Since the mass range of interest in thrustor design lies between 10 4 and 10 6 
atomic mass units (ref. 2), aggregates containing 10 3 and 10^ atoms for each of 
of the alkalies, lithium, sodium, potassium, cesium, and rubidium, are con- 
sidered. 
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CALCULATION OF IONIZATION TIME 

The electrons in a metallic colloid move under an average potential that is 
lower than the vacuum potential. At the surface, a potential 11 jump 11 to the 
vacuum level restrains the electrons from leaving the colloid. The probability 
of finding an electron at an energy greater than the vacuum level is exceedingly 
small at ordinary temperatures, as can be seen from the Fermi distribution given 
by equation (l). 

An external electric field Impressed on the colloid creates a potential 
barrier through which electrons can tunnel. The potential at the surface of the 
colloid is at the vacuum level and decreases linearly with a slope proportional 
to the field strength. The number of electrons n passing through this barrier 
per unit time depends on the number striking the barrier per second and on their 
transmission probability through the barrier. The reciprocal of this number is 
the average time required for one electron to tunnel through and is defined as 
the charging time t: 


T 



( 2 ) 


If the field is in the x-direction and dn is the number of electrons within 
the colloid having the velocity v x and a transmission probability D(Ex), then 
it follows that 


<in = Sf (E)v x D(E x )dp x dp y dp z (3) 

where 2/h^ is the density of levels (two spin states per volume h^ of phase 
space, h being Planck* s constant); S is the effective surface area normal to 
the x-direction; and p x , py, and p z are the momentum components. The trans- 
mission probability is only a function of the energy in the x-direction. 


The total number of electrons that will penetrate the barrier is then found 
by integrating over all possible momentums of the y- and z-components and over 
all positive values of the x-component: 
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When the last two integrations are carried out, there remains 
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Before this integration can be completed, it is necessary to find the trans- 
mission factor D(E X ) as a function of energy. Because the potential is only a 
function of the x- coordinate, the wave equation is reduced to a one -dimensional 
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problem, and hence the W.K.33. (Wentzel, Kramers, and Brillouin) approximation 
(ref. 4) is applicable. Accordingly, 



•where the integration is carried out over the barrier. To simplify the calcula- 
tion, V(x) is approximated by a triangular barrier depicted in the following 
sketch: 



x 


(a) 

The limit of the harrier at x = a is seen from the sketch to be a function 
of both the applied electric field F and the energy Ex of the impinging elec- 
tron. Carrying out the integration gives 

D(E i> - - ^^heF 172 (E ® ' E * )3/2 (7) 

The ionization time is calculated with equation (7) substituted for 
D ( E ) in equat ion ( 5 ) : 



( 8 ) 

Since this integration cannot be carried out in closed form, an approxima- 
tion will be sought. Most of the contribution to the integral comes from below 
the Fermi level; therefore, the first simplification is to replace the upper 
limit by E Q and to let 
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Next, by expansion of (Eg - E x ) ‘ in a Taylor series about E x = E Q and 
by retention of only the first two terms, an integrable approximation is ob- 
tained: 
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where cp, the work of function of the material, is equal to Eg - E Q . 

Because the transmission factor rapidly converges to zero, the lower limit 
in equation (8) can be replaced by -<». With the substitution t) = (E q - E x )/cp, 
the reciprocal charging time finally takes the form 
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Fowler and Nordheim 1 s equation (ref. 5) for the current density emitted 
from a metal in a strong electric field shows the same functional dependence on 
the field strength as equation (9) although the method of derivation was not the 
same. 


The effective surface area S in the field direction is taken to be the 
cross-sectional area of a spherical droplet of volume V: 


S = 



it 



( 10 ) 


Where N is the number of atoms in the colloid and p a is their number density, 
taken to be equal to the bulk density. Figure 1 shows a plot of charging time 
as a function of field strength for alkali-metal colloids containing 10 ^ and 10 ^ 
atoms. Table I contains the information used in calculating the curves shown in 
the figures. 
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IONIZATION PROBABILITY AND UTILIZATION EFFICIENCY 


In order to produce the strong electric fields needed for ionization, it is 
necessary to use small geometries such as needle points or thin vires. An added 
effect of the field is produced by the neutral colloids being dravn into the 
regions of high field strength because of polarization forces. The polarization 
force afforded by the field of a cylinder is a longer range effect (~r~ 3 ) than 
the field of a point (~r“5). p or this reason, the folloving geometry, consist- 
ing of three coaxial cylinders, is considered: 



The central wire and the outermost cylinder are maintained at high potentials 
vith respect to the grounded cylindrical grid. This arrangement forms a "poten- 
tial veil," as shovn in the sketch. 

It is assumed that a veil collimated noninteracting beam of colloidal par- 
ticles is introduced axially into one end of the tube. Because of the polariza- 
tion forces, the particles are gradually dravn into the high-field region sur- 
rounding the central vire, charged, and repelled into the lover potential region. 
The charged particles oscillate in the potential veil until they drift out of 
the tube. Once they are out, they can be introduced into a suitable thrustor. 

The classical expression for the polarizability a of a spherical con- 
ductor of radius R is 


a = 4jte Q R 3 (ll) 

This becomes for spherical colloids 

a = 3e 0 ~ (12) 

° Pa 

where N and p , as before, are the number of particles and the density of the 
droplet, respectively. 
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The potential energy of such a colloid in an electric field is then 


U = 



3 

2 



( 13 ) 


The field about the wire is radially symmetric and can be written in terms 
of the wire radius r a and the surface field strength F a : 


F 



(14) 


Since the beam is assumed to be well collimated, the angular component of 
the particle velocity is ignored, and the radial component is obtained from the 
equation of motion 


dr 

dt 



(15) 


where m* is the atomic mass as given in table I (a) and F^ is given by equa- 
tion (14) at r - r j_, the initial distance from the origin. 


The z-component is not affected by the field and is conserved through the 
tube even after ionization if no particle interactions occur. 


The probability P that the colloid will emit an electron as it moves into 
the high-field region is given by 



where dt is the differential residence time from equation (15). In terms of 
the field strength, the integral in the exponent becomes 



The charging time t in equation (17) is a function of an enhanced field 
strength, which has an average value of three-halves of the external field. 

This arises because of the distortion caused by the polarization of the colloid. 

Initially, the colloids can enter the field at all possible values of Fj_; 
that is, at any distance r^ from the anode. The integrand, however, is a very 



sensitive function of the field strength and rapidly vanishes away from the 
anode. For this reason, the initial field strength is neglected. This is 
equivalent to carrying out the calculation for a colloid starting at infinity. 
Colloids nearer the anode will, of course, also become ionized under the same 
conditions. 

If the substitution 


x 



9heF r 
a a 


(18) 


is made, the integration takes the form of the well-known exponential integral 
-Ei(-x a ), which is tabulated in almost all mathematical handbooks (e.g., ref. 6). 
In the expression 


-Ei(-x a ) 



(19) 


x a is obtained from equation (18 ) with r = r a . With these substitutions and 
definitions, equation (17) becomes 



The ionization probability (eq.. (16)) is very sensitive to small changes in 



These same changes, however, are reflected in equation (20) as a slight 


modification of the field strength (for a given radius); that is, the same order 
of magnitude field strength is needed to produce partial ionization of the beam 
as is used to achieve almost total ionization. In view of this, a condition 


will arbitrarily be chosen by letting 



which yields a probability of 


0.982. The utilization efficiency at this degree of ionization of the beam 
should then be greater than 98 percent provided that the neutrals do not leave 
the tube before they come into the ionization region. Solving equation (20) for 
r a with a value of 0.982 for the integral then gives the desired relation for 
the wire size as a function of the necessary field strength: 
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The. potential V a between the wire and the grounded grid cylinder of radius r Q 
is given by 


T a * F a r a ln < 2 2) 

The field strength can be eliminated by combining equations (21) and (22) 
to obtain the voltage requirement as a function of wire radius. These are two 
of the more important design parameters plotted in figure 2 for grid radii of 
1CT 3 and 10" 2 meter. 


To complete the discussion of utilization efficiency, it now remains to 
calculate the necessary cylinder length to allow each colloid to reach the ion- 
ization zone before leaving the tube. 


The fraction of the colloidal beam having a velocity between v z and 
v_ 4 - dv„ is given by 


f O )dv 
v z 7 z 

Integrating from zero to v z 
less than v z : 
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(23) 


then gives the total fraction having velocities 


f = 1 - (l + v 2 ) exp (-v 2 ) 
z 


(24) 


where 


v 



(25) 


The time it takes a particle starting at the beam edge r-j_ to reach the 
anode (r a ~ o) should be made equal to its transit time through the tube travel- 
ing at a velocity v z . Hence,, by use of equation (15),, 
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from which 


L = 



l/2 vr| In 


V 

a 



(27) 


Table II gives the value of L at the atmospheric melting point of each of 
the alkalies at the following parameter values: 


N 

xH 

O 

H 

V 

1 

r Q = r ± , m 

10~ 2 

V m 

LO 
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O 
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10 5 


If this value Is represented by L*, equation (27) can then be rewritten in terms 
of its more important parameters as 

r 2 T * 

L = vp^L 


where pq is a nondimens ional factor equal in magnitude to the radius of the 
beam expressed in millimeters, and v, of course, gives the percent of beam f v 
that will be Ionized as determined from equation (24). z 

When the values of r D , r a , and V a were chosen, consideration was given 
to figure 2, which shows these values as occurring at realistic operating points. 
Except for variations in v and p-j_, equation (27) must be used as the basis for 
computing the tube length for other values of these parameters. 


CHARGE AND CURRENT DENSITIES 

The purpose of maintaining the outer cylinder at a high potential is to 
provide a sufficiently strong decelerating force on the charged colloids to keep 
them confined within the tube. In order to accomplish this, the potential must 
be at least as high as the potential on the emitting wire. 

The charged particles are formed in the high-potential region around the 
wire. From figure 2, it is evident that, for practical emitter radii and oper- 
ating potentials, each colloid would have an energy of the order of 10^ electron 
volts if singly charged particles are assumed. The maximum density of neutral 
particles entering the tube is estimated to be of the order of 10^-5 per cubic 
centimeter. The interaction energy of two charged colloids at this density 
would then be of the order of 0.01 electron volt. The description of the motions 
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of the charged colloids 'within the tube then coincides 'with that of an ideal gas 
in 'which the particles move independently of each other hut achieve thermodynamic 
equilibrium by virtue of the randomizing effect of the weak interactions be- 
tween particles. 


With the Boltzmann factor, the relative density of particles within the 
tube can be found. The potentials in the inter cylindrical space, for the case 
where the outer cylinder is at the same potential as the wire, are given as 


v(r) 
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In its most general form, the density can be written as 

H 


P = c 


/ • / exp fii) dp i 


dP 


3N 


(29) 


where H(p^ . . . q^ . . . q.^) is the Hamiltonian function of the conjug- 

ate momentums and coordinates p and q, and c can be considered to be a 
normalizing constant. Since the potentials depend only on the coordinates, the 
integrations over the momentums contribute a constant factor and the density is 
a function of the coordinates. In particular, because of equation (28), the 
density is a function of the radial coordinate alone: 


p(r) = p(r Q ) exp 


eV(r) 

kT 


(30) 


where p(r 0 ) is a constant equal to the maximum density that occurs at V = 0, 
that is, at r = r Q . In terms of the potentials of the two regions, as given in 
equation (28), the normalized density is then 


p 





< r < r 0 




(31) 
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The quantity eV a /kT needed in equations (3l) requires a knowledge of the 
ion gas temperature. Since the initial conditions are approximately known, the 
effective temperature can he calculated. The average energy per emitted particle 
is approximately eV a , while from the ideal gas equation, the average energy is 
3 

— kT. Equating these terms gives the desired relation 



kT 


3 

2 


(32) 


A field- emission ionizer, such as the one described, can never operate under 
space -charge -limited conditions, for this implies that the field at the emitter 
vanishes. It is not clear what the exact saturation level is, but it is obvious 
that the field does not have to be suppressed very greatly to cause a notable 
increase in the charging lifetime and thus end field emission. If it is assumed 
that an anode field-strength suppression of 10"^ F a is the allowable maximum, 
an estimated upper limit to p(r Q ) (appendix B) is 


p(r c ) 



(33) 


The average flux density, out of the tube end, with an assumed mean 


drift velocity v , is then equal to 
z 


T = — 

J E 2 
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The mean velocity of a molecular beam in the mass range being considered 
(Em* ~ 10“^ kg), at a temperature of the order of 1000° K, is approximately 
7 meters per second, as computed from 


v 


z 


3 /2rtkT \ 1//2 

4 \Nm* / 


(36) 


V 


a 


The value of calculated at a likely design point, for instance, 

cr ^ * cr 'z 

10° volts, r a = 10”° meter, and r Q = 10" meter is therefore of the order 
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of 10 1 ® per square meter per second. This figure, if multiplied by the elec- 
tronic charge e gives a maximum current density of 0.144 ampere per square 
meter from this device in thermal equilibrium; that is, if no field is used to 
draw the charged colloids out of the tube. This figure also represents the 
necessary input of neutral colloids from the particle generator, in order to 
achieve maximum charging efficiency More than this would cause saturation and 
a resulting loss in ionization efficiency. The use of even a moderate electric 
field applied axially along the cylinder, however, can greatly increase the cur- 
rent density. 


CONCLUDING REMARKS 

The field emission method of producing positively charged colloids is, on 
the basis of net power per charge, one of the most efficient ways of charging. 

The propellant -ut ilizat ion efficiency, as pointed out in the text, could be 
made to approach unity in a cylinder of reasonable length. This, of course, is 
contingent on how well the colloid generator can form aggregates from the mona- 
tomic propellant. 

At present, the necessary relation between the applied potentials and the 
emitter radii required for this charging method is somewhat difficult to 
achieve, but it is likely that in the future monocrystalline filaments and/or 
high-voltage sources will be available. With the advent of such thin filaments, 
it is possible that the required emitter potential will be reduced enough so 
that a magnetic field can be used to turn the charged colloids, and thus the 
intermediate grid could be eliminated. This is desirable since the grid is 
located in the region of greatest charge density and is subjected to severe im- 
pingement. 


Lewis Research Center 

National Aeronautics and Space Administration 
Cleveland, Ohio, August 13, 1963 
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APPENDIX A 


a 

c 


D 

E 


E. 


B 


V E z 



N 

No 


SYMBOLS 

potential barrier end point 
constant 

transmission probability- 
energy level or state 
vacuum potential 
Fermi energy level 

translational energy in x~, y-, and z -directions, respectively 

electronic charge 

electric field strength 

electric field strength at anode 

electric field strength at r^_ 

Fermi distribution function; velocity distribution function 

fraction of particles having velocities up to v_ 

z 

Hamiltonian function 
Planck' s constant 
average flux density 
Boltzmann' s constant 
cylinder length 

specific cylinder length, and table following it defined by 
eti. (28) 

electronic mass 

atomic mass 

total number of constituent atoms in colloid 
Avogadro's number 


14 



n 

On 


P 

P 

W P z 

q. 

R 


r 


r i 

r o 

S 

T 

t 

U 

V 


Vb 

v x> v y> v z 
x 


x a 


*>y, z 


r, 9, z 


member of electrons passing through potential harrier per sec 

differential number of electrons passing through potential barrier 
per sec 

ionization probability 
momentum 

momentum components in x-, y-, and z-directions, respectively 
coordinate 

effective colloid radius 
radius 

anode radius 

outer cylinder radius 

initial radial distance from anode 

grid radius 

effective surface area presented to tunneling electrons 

temperature 

residence time 

potential energy 

potential function 

voltage at anode 

voltage at outer cylinder 

velocity components in x-, y-, and z-directions, respectively 
average velocity in z-direction 
dummy variable in eq_. (20) 
defined by eq.. (19) at r = r & 

Cartesian coordinates 
cylindrical coordinates 
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a 



a 

T 


9 


polarizability 
vacuum dielectric constant 
dummy variable 
nondimensional beam radius 
variable defined by eq. (27) 
density function 

number density of atoms in colloid 
density of water at 4° C 
surface variable 
charging time 
work function 


Subscripts: 


i initial 

1 . . . 3N degrees of freedom 
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APPENDIX B 


SUPPRESSION OF ANODE FIELD STRENGTH 


The field strength at a distance r due to a cylindrically symmetric cloud 
charge of density ep(r) and a central wire conductor having a surface charge 
density of e F(r_) is 

O 0, 



(B1 


The -wire is held at a constant potential V a with respect to the grounded 
grid at r Q . Hence, this voltage is given alternately "by equations (22) and 
(B2): 


v a = 


-/ 


F(r)dr 


(B2) 


If the appropriate expression for p(r) from equation (31) is used and the 
indicated integrations are carried out, the effect of the space charge on the 
wire field strength is found to he 


F *. - F(rJ = 


e p( r 0 ) r o 


€.r In _£ 
o a y» 


2 ln(i- 0 /r a ) 


+ 2 


Terms involving powers of r a are ignored since r D » r a . 
If a field-strength suppression factor of 

F a - F ( r a) = 10 ” 2 F a 


is assumed, the maximum density becomes 




P ( r o ) = 


I 


2 ln(r c /r a ) 


"l ‘ i 

+ 2 


er r 


(B3) 


(B4) 


(B5) 


in which equation (22) was again used. 
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TABLE II 


CYLINDER LENGTH AT 


MELTING POINT 


Element 

Melting point 
at 760 mm Hg, 
°K 

Specific 

cylinder 

length, 

L* , 
m 

Lithium 

459.2 

3.25X10 -1 

Sodium 

370. 7 

2. 19X10“ 1 

Potassium 

335.5 

1.48X10 -1 

Cesium 

301.7 

3.54X10" 2 

Rubidium 

311.7 

4. 10X10" 2 









Anode radius, r a/ m 

(b) Atom density, 10 4 atoms per colloid; grid radius, 10“ 2 meter. 
Figure 2. - Continued. Anode voltage as function of anode radius 




Anode radius, r a , m 

(c) Atom density, 10^ atoms per colloid; grid radius, 10“^ meter. 
Figure 2. - Continued. Anode voltage as function of anode radius 
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(d) Atom density, 10^ atoms per colloid; grid radius, 10"^ meter 


Figure 2. - Concluded. Anode voltage as function of anode radius 



